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Determining the properties of the dark components of the universe remains one of the outstand- 
ing challenges in cosmology. We explore how upcoming CMB anisotropy measurements, galaxy 
power spectrum data, and supernova (SN) distance measurements can observationally constrain 
their gravitational properties with minimal assumptions on the theoretical side. SN observations 
currently suggest the existence of dark matter with an exotic equation of state p/p < —1/3 that 
. . . accelerates the expansion of the universe. When combined with CMB anisotropy measurements, SN 

OO , or galaxy survey data can in principle determine the equation of state and density of this component 

separately, regardless of their value, as long as the universe is spatially flat. Combining these pairs 
creates a sharp consistency check. If p/p ^ —1/2, then the clustering behavior (sound speed) of the 
dark component can be determined so as to test the scalar-field "quintessence" hypothesis. If the 
' exotic matter turns out instead to be simply a cosmological constant (p/p = —1), the combination 

^ ' of CMB and galaxy survey data should provide a significant detection of the remaining dark matter, 

the neutrino background radiation (NBR). The gross effect of its density or temperature on the ex- 
' pansion rate is ill-constrained as it is can be mimicked by a change in the matter density. However, 

\ anisotropies of the NBR break this degeneracy and should be detectable by upcoming experiments. 
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I. INTRODUCTION 



> 

' The nature of the dark matter remains one of the greatest outstanding puzzles in cosmology. The difficulty of the 
problem is compounded by the fact that the dark matter may be composed of multiple components. Despite this, we 
may be on the verge of an observational solution. The cosmic microwave background (CMB) contains information 
about the dark components present in the early universe, specifically the ratio of non-relativistic or cold dark matter 
(CDM) to relativistic species such as the neutrino background radiation (NBR) and the ratio of the baryonic dark 
matter to the CMB itself. Upcoming high precision measurements of the CMB, notably by the MAP |jl[| and Planck 
7*^ 1^ satellites, should determine these ratios to the percent level In contrast, observations of high-redshift objects 
O such as Type la supernovae (SN) probe dark components important in the local universe. Indeed preliminary results 
> ■ suggest the presence of an additional dark component that accelerates the expansion . The clustering properties 
^ ' of galaxies link the CMB and the local universe through their dependence on both the initial perturbations visible in 
!• i: the CMB and the time- integrated history of structure formation between last scattering and the present. The galaxy 
. 5^ ' power spectrum will be precisely measured by ongoing redshift surveys such as the 2dF |^ and the Sloan Digital Sky 
>^ . Survey (SDSS) §. 

5h ] The promise of observationally determining the properties of the dark components lies in combining these data sets. 
. - - ' Aside from the obvious difference in redshift windows, the various data sets individually suffer from the fact that their 
observables depend degenerately on several aspects of the cosmology. Combined, they break each other's degeneracies. 
With three data sets, consistency tests become possible. These tests are valuable for investigating systematic errors 
in the data sets and could potentially indicate that the current cosmological framework is inadequate to describe the 
universe. 

In this paper, we investigate how the combination of CMB anisotropy measurements, galaxy survey data, and SN 
luminosity distance determinations can be used to determine the parameters of the dark components. For this purpose, 
we employ the generalized dark matter (GDM) parameterization scheme introduced in ^J. This parameterization 
encapsulates the observable properties of the dark components in a background equation of state, its density today, 
a sound speed, and an anisotropic stress or "viscosity parameter" . We begin by examining how well the equation of 
state and density today can be determined from observations assuming a flat universe. Combining CMB data with 
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either galaxy surveys or SN observations will provide tight constraints on the equation of state and the density of the 
exotic component even if the sound speed or viscosity must be simultaneously determined. The combination of these 
pairs will thus provide a sharp consistency test. Galaxy survey information assists in these measurements indirectly 
by freeing CMB determinations from parameter degeneracies. Its power is revealed only upon a full joint analysis and 
cannot be assessed by direct examination of individual parameter variations (c.f. [^). The fundamental assumption 
is that the galaxy and matter power spectra are proportional on large scales where the fluctuations are still linear. 

As the cosmological constant has a well-defined equation of state p/ p = —1, these cosmological measurements will 
test for its presence. A cosmological constant is also special because its density remains smooth throughout the 
gravitational instability process. If the exotic component proves not to be a cosmological constant, then its clustering 
properties become important. These properties are encapsulated in the sound speed. The simplest models for this 
component involve a (slowly- rolling) scalar-field "quintessence" pO|-p^ which has the interesting property of having 
the sound speed in its rest frame equal to the speed of light ]14|,|8| . By measuring the sound speed, one tests the 
scalar-field hypothesis. As long as p/ p ^ —1/2 in the exotic component, the combination of CMB experiments and 
galaxy surveys can provide interesting constraints on the sound speed. Finally, we show that if the exotic component 
turns out to be simply a cosmological constant, then the properties of the remaining dark matter, the NBR, can be 
determined from combining CMB and galaxy survey data. In particular, anisotropics in the NBR, as modeled by the 
viscosity parameter of GDM, are measurable and provide a way to break the matter-radiation density degeneracy in 
the CMB. 

The outline of the paper is as follows. We review the phenomenology of the parameterized dark matter model of |q| 



n] and the Fisher matrix technique for parameter estimation in §[II. In §IV, we determine how well the equation 



of state of the dark component may be isolated from its density. We discuss measurements of the sound speed and 
propose a test of the scalar-field quintessence hypothesis in §0. In § VI , we address the detectability of the NBR and 



its anisotropics. We summarize our conclusions in §VII 



II. DARK MATTER PHENOMENOLOGY 



In this section, we review the observable properties of the dark components following the phenomenological treat- 
ment of We shall see that the properties of the dark sector affect CMB anisotropies, structure formation, and 
high-redshift observations in complementary ways. 

We define the dark sector to include all components of matter that interact with ordinary matter (baryons and 
photons) only gravitationally. Thus, the observable properties of the dark sector are specified completely once the 
full stress-energy tensor is known. Here we consider the dark sector to be composed of background radiation from 
3 species of essentially massless ^ O.leV) neutrinos, a cold dark matter (CDM) component, and an unknown 
component of generalized dark matter (GDM) The CDM is required to explain the dynamical measures of the 
dark matter associated with galaxies and clusters. The GDM is required to be smooth on small scales to avoid these 
constraints [ p^ljTsI] . We further assume these forms of dark matter do not interact at the redshifts of interest. 

Since each non-interacting species is covariantly conserved, the ten degrees of freedom of the symmetric stress-energy 
tensor of the GDM are reduced to six. We take these as the six components of the symmetric 3x3 stress tensor. Two 
stresses generate vorticity and two generate gravity waves; we will not consider these further but note that they may 
play a significant role in CMB anisotropy formation in so-called "active" models for structure formation This 
leaves two stresses: the isotropic component (pressure) and an anisotropic component ("viscosity"). The properties 
of these two stresses must be parameterized. We begin by discussing their effect on the background expansion and 
then examine their role in the gravitational instability of perturbations. 



A. Background Effects 

Isotropy demands that to lowest order the GDM stress tensor has only an isotropic (pressure) component. The 
GDM properties to lowest order therefore depend only on the equation of state Wg = Pgj pg. For example, energy 
conservation requires that the evolution of the GDM density follows 

dln(l+Vi =-^(^ + "^)- 
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FIG. 1. Distance measure degeneracies: contours of constant (a) luminosity distance to z — 0.5 (HodL) and (b) angular 
diameter distance to the last scattering surface (HocIa)- ^mh? has been held fixed in the former under the assumption that 
the CMB acoustic peak morphology will measure it independently. With this assumption, the two distance measures provide 
complementary information. 



For simplicity, we consider models where Wg is independent of the redshift z since the properties of a slowly-varying 
Wg can be modeled over the relevant redshifts with an appropriately weighted average [^|T|. Combined with the 
assumption of zero spatial curvature, the expansion rate or Hubble parameter becomes 

^{ng,wg,n,r.inr;z) = [n„,{i + zf + ng{i + zf(^+^^) + nr{i + zfY'\ (2) 

no 

where Hq — 100 /i km s^^ Mpc^^. Here Vl„i = 1 — Q.g — Vlr and accounts for both the CDM and baryonic ("matter") 
components. Likewise £7^ accounts for the photon and neutrino ("radiation") components. Furthermore, quantities 
that depend on the redshift behavior of the expansion rate such as the deceleration 

^ = H d(l + z)-i+^ 

depend on the same parameters. In particular, a component with Wg < —1/3 is able to drive q negative and cause an 
acceleration. 

Any cosmological observable that is simply a function of the expansion rate of the universe will only be sensitive to 
the background properties of the GDM through {wg, fig). For example, SN probe the luminosity distance (see Fig. |^a) 

HodL{wg,ng;z)^{l + z) f dzHo/H; (4) 

Jo 

this function is independent of the matter-radiation ratio because the observations are at sufficiently low redshift. 
Current SN luminosity distance data suggest the presence of an accelerating component with Wg < —1/3 (see Fig. H) 
but may be dominated by unknown sources of systematic errors. 

Likewise the acoustic peaks in the CMB probe the angular diameter distance to the redshift of last scattering zis 
(see Fig. |b) 

HodA{wg,ng,n„i/Qr;zis) = / dzHo/H, (5) 

Jo 

through its ratio with the sound horizon at last scattering HoScMB{^-y/^b,^m/^r)- Current CMB detections alone 
do not place significant on Wg and ng [H. 

The energy density of the CMB is fixed through the measurement of its temperature Tfiras = 2.728 ±0.004K (95% 
CL) H such that 



^ 2.4815 X 10-^ (^^\ 

V ^FIRAS / 



(6) 
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FIG. 2. Current SN data. The 65%, 95% and 99% CL intervals in the u 
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- fig plane for the current data assuming only 
statistical errors. Constraints include 6 high redshift SN from the Supernova Cosmology Project Q and 10 from the High-z 
Supernova Search We use 26 low-2 calibrating SN with B —V < 0.2 obtained by the Calan/Tololo group [^. The analysis 
follows |pi[, but note that systematic errors may dominate in the current data sets. 



With this constraint, the dependence of the sound horizon on the photon-to-baryon ratio reduces to a dependence on 
flbh'^. Likewise, the total radiation energy density is given by 
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(7) 



Thus under the usual assumptions for the number of neutrino species {N^ = 3) and their thermal history {T^ = 
101^16/1.401), the dependence of cZa and scmb on flm/flr becomes a dependence on Qmh^- We relax these assumptions 
in § [V^ to test the properties of the NBR. 



B. Structure Formation 



To probe the remaining properties of the GDM, one must consider its effects on the gravitational instability of 
perturbations. Unless Wg = —1 (the cosmological constant case fig = JIa), the GDM participates in the gravitational 
instability process. The cosmological constant case is special since the relativistic momentum density, which is 
proportional to 1 + Wg, goes to zero. 

Since the perturbations need only be statistically isotropic, the GDM in general requires two parameters to describe 
fluctuations in its stress tensor. These can be chosen to be the sound speed in the rest frame of the GDM CeS, where 
c^jj = Spg/Spg (in units of where c = 1), which relates the pressure fluctuation to the density perturbation, and a 
"viscosity" parameter Cvis, which relates velocity and metric shear to the anisotropic stress. See Q for their precise 
covariant definition. Positive values of c^g imply that density fluctuations are stabilized by pressure support at the 
effective sound horizon Scff = / Ccfr(l + 2)<ii. Likewise, positive values of c^jg imply that resistance to shearing stabilizes 
the fluctuation at Svis = /cvis(l + z)dt. These definitions assume that Ccs and Cvis, respectively, are slowly varying. 
We call the greater of ScS and Svis the GDM sound horizon sgdm- 

Modes smaller than sgdm are stabilized by stress support. If the GDM also dominates the expansion rate then the 
growth of structure will slow below this scale. If Wg < 0, GDM domination occurs at approximately 



4 




FIG. 3. Sound speed effects for = — 1/3. A finite sound speed Cefi stabilizes perturbations leading to features in the 
(a) galaxy power spectrum and (b) CMB anisotropy spectrum. Note that in both cases the effects change most rapidly with 
Ccff between and ^JlJ&. The power spectra have been normalized to small scales to bring out the degeneracies and the 
importance of large-scale information. The model here and throughout has Q,m = 0.35, h = 0.65, Q,f,h? = 0.02, r — 0.05, n = 1, 
and T/S = 0. 



Thus we expect a feature in the matter power spectrum between i?oSGDM(wg, Vlg) ai z = Zg and z = 0. Since below 
the sound horizon, the effect of GDM is to slow the growth of structure independent of scale, the determination of Coff 
from measurements of the galaxy power spectrum depends crucially upon having data across this range of scales. In 
Fig. we show the effect of varying Coff on the power spectrum for Wg = —1/3. The models have been normalized 
to small scales to bring out the scale-independence of the small-scale suppression. Given that the normalization is 
uncertain because of the unknown proportionality constant between the mass and galaxy power spectrum usually 
defined as 6^ where b is the "bias", the only direct information from galaxy surveys on CcS comes from large scales. 
Notice the most rapid change with Coff occurs for CcS < \/l/6. 

The location of the GDM sound horizon also affects CMB anisotropics. By halting the growth of structure, the 
GDM causes gravitational potential wells to decay below the sound horizon sqdm &t Zg. A changing gravitational 
potential imprints fluctuations on the CMB via differential gravitational redshifts whose sum is called the integrated 
Sachs- Wolfe (ISW) effect. However, if the sound horizon is much smaller than the particle horizon, the photons will 
traverse many wavelengths of the fluctuation as the potential decays. The cancellation of red and blue shifts destroys 
the effect. Thus the largest effect arises when Coff ~ 1 but varies strongly with Coff only as it becomes less than 
Coff ^ -\/l/6. Unfortunately, subtle differences in this large-angle temperature signal around c^s — 1 will be difficult 
to pin down given cosmic variance. 

Because of the strong Ceff-dependence of features in the CMB and galaxy survey power spectra for Ceff % -^/l/G, 
tight lower limits can be placed on Ccff from the data even though models around Ccg = 1 cannot be distinguished. 
Furthermore, the amplitude of the features decreases sharply as —1 since even clustering above the sound 

horizon vanishes in this limit. Constraints on the sound speed will thus only be possible if the equation of state of 
the GDM differs significantly from a cosmological constant. 

Finally, the viscous term also causes the growth of structure to halt. However, it has an additional effect that 
makes it unique since anisotropic stresses enter directly in the Poisson equation that defines their relation to the 
gravitational potentials . In Fig. ^, we show the effect of replacing the neutrinos with GDM of Cvis — ^^1/3 and 
in the CMB. The former component models the neutrinos accurately and the latter shows that the anisotropic stress 
of the dark matter produces potentially observable effects in the CMB. We will exploit this effect to propose a means 



of detecting the anisotropics in the neutrino background radiation in § VI . 



III. PARAMETER ESTIMATION 



Projections for how well various data sets can measure cosmological parameters depend crucially on the extent of 
the parameter space considered as well as on the location in this space (or "fiducial model" ) around which we quote 
our errors. Even though the latter uncertainty will be eliminated once the best fit model is found from the actual data, 
the former problem will remain. The extended GDM parameter space thus allows the data the freedom to choose the 
best values to describe properties of the dark sector. Even if the true model turns out to contain only conventional 
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dark matter, it allows us to say with what confidence we can make this statement, i.e. that the dark components are 
in fact the neutrino background radiation, a cosmological constant, and cold dark matter. If these options are ruled 
out, then we will have discovered a new form of matter. 

We adopt a 10-dimensional parameterization of cosmology that includes the present density of the GDM fig, a 
time-independent equation of state Wg, the matter density flmh^, the baryon density Q.hh'^, the reionization optical 
depth T, the tilt n, the tensor-scalar ratio T/S^ the normalization A, and the linear bias b. Both c^s and Cvis affect 
the clustering scale and are largely degenerate. In §[V and we take Coff as a proxy for both, whereas in §VI, we 
take Cvis since we are interested in the anisotropy itself. This parameter space does not include models with non-zero 
spatial curvature or with massive neutrinos. We take the fiducial model to have r2„j — 1 — ^Ig — 0.35, h — 0.65, 
fife/i^ = 0.02, r = 0.05, n = 1, and T/S = 0; we will use fiducial models with different values of Wg and Ceg to explore 
how the results depend on these parameters. For how the fiducial choices in the standard parameters affect the errors, 
see|,||. 

To estimate errors on the cosmological parameters, we employ the Fisher matrix formalism (see for a general 
review). The Fisher matrix is essentially an expansion of the log-likelihood function around its maximum in parameter 
space. It codifies the optimal errors on each parameter for a given experiment assuming a quadratic approximation for 
this function. Note that the Fisher matrix provides accurate estimates of error contours only when they encompass 
small variations in parameter space where the expansion is valid. This is an important caveat to which we will return 

Fisher matrix errors thus employ derivatives of the cosmological observables with respect to parameters. To calculate 
derivatives of the CMB and galaxy power spectra, we employ the hierarchical Boltzmann code of ||2^,^. Special 
care must be taken in evaluating these derivatives since numerical noise in the calculation can artificially break any 
parameter degeneracies that exist. General techniques such as the taking of two-sided derivatives and the associated 
step sizes for standard parameters are described in ||2^ . For the GDM parameters, we estimate the derivatives by 
finite differences with the step sizes Afig — ±0.05(1 — Ug), Awg — ±0.01, Ac^fj = ±0.1clg and Ac^j^ = ±0.1c^jg. 

The benefit of a hierarchy treatment is that unlike the integral treatment of CMBfast no interpolation is 
necessary and the code may be made arbitrarily accurate by adjusting the sampling in Fourier space. Thus the 
numerical noise problems identified by p3t can be addressed and in principle eliminated. However, computational 
speed generally requires a compromise involving smoothing the calculated CMB power spectrum p6t . We have tested 
our results against CMBfast version 2.3.2 in the context of flat A-models where CMBfast is most accurate and obtained 
~ 5% agreement in parameter estimation with the hierarchy code. The test also involved two independent pipelines for 
taking model calculations through to parameter estimations. These comparisons were done with 1500 Fourier modes 
out to A: = 2.5A:damp, where fcdamp is the CMB damping scale (see |Q, Eq. [17]) with Savitzky-Golay smoothing |2|] 
of the resulting CMB angular power spectrum. We adopt the same techniques for parameter estimation in the GDM 
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context .[| 

In addition to the extent of the parameter space and the location of the fiducial model in parameter space, Fisher 
errors of course depend on the sensitivity of the given experiment. For the CMB data sets, we take the specifications 
of the MAP and Planck experiments given in . We quote results with and without polarization information. Since 
the polarization signal may be dominated by foregrounds and systematic errors, the purely statistical Fisher matrix 
errors may be underestimates. For galaxy surveys, we take the Bright Red Galaxy sample of SDSS; the specifications 
and their translation into the Fisher matrix formalism is given in pO[ . We further take the linear power spectrum 
for parameter estimation. Because non-linear effects and galaxy formation issues complicate the interpretation of the 
observed power spectrum on small scales, we employ information only from wavenumbers less than fcmax = 0.2/iMpc~^ 
and show how the results change as we go to a more conservative fcmax = 0.1/iMpc~^. This roughly brackets the 
regime where non-linear effects begin to play a role as shown by simulations For SN, we assume that on a 

timescale comparable to the MAP mission a total of 200 supernovae will be found with individual magnitude errors 
of 0.3 and a redshift distribution with a mean of 2 = 0.65 and a gaussian width of Az — 0.3 Wh - 

As each data set is independent to excellent approximation, the combined likelihood function is the product of 
the CMB, SN and galaxy survey likelihood functions. Thus the combined Fisher matrix is simply the sum of the 
individual Fisher matrices. 



IV. MEASURING THE EQUATION OF STATE 

Current SN luminosity distance measures suggest that the GDM may have an exotic equation of state Wg ^ ~l/3 
that accelerates the expansion (see Fig. |2]and If these preliminary indications are borne out by future studies, 

one would like to pin down the equation of state of the GDM and also construct consistency tests to verify this 
explanation of the SN data. Unfortunately, no one data set can isolate the equation of state on its own. The 
CMB has one measure of Wg from the angular diameter distance to the last scattering surface (see Eq. j^), but this is 
degenerate with f2g (see Fig. |^), even assuming zero curvature and that flmh^ has been measured from the morphology 
of the acoustic peaks. Leverage on Wg and Slg comes only through the cosmic-variance-limited ISW effect at large 
angles and through the effects of gravitational lensing on very small scales. The latter occurs because changing fig 
affects the present-day normalization of the matter power spectrum and thereby changes the amount of lensing. The 
effect is small, however, and much less powerful for breaking the degeneracy than the methods of the next paragraph; 
we therefore neglect lensing in the CMB. Likewise, galaxy surveys alone give leverage only through the combination 
that defines the GDM sound horizon. For SN measurements that span only a short range in redshift, there is an 
analogous degeneracy between Wg and fig in the luminosity distance (see Eq. Q). 

However, combining two of these measurements isolates Wg and the third can be used as a consistency check. That 
the CMB angular diameter distance and SN luminosity distance measures break each others degeneracies is obvious 
from comparing panels (a) and (b) in Fig. |l]. We show in Fig. ^ the error ellipses (68% CL0) in the Qg — Wg plane. 

Due to the large ISW effect in the Wg — —1/6 model of Fig. ||a, MAP alone will provide reasonable constraints 
on the two parameters. In this case, SN measurements will provide a strong consistency check on CMB measure- 
ments. However, as one approaches Wg = —1, the ISW effect decreases and the CMB requires the assistance of SN 
measurements to break the degeneracy. We list the la errors as a function of Wg in Tab. ^ 

The combination of CMB and galaxy survey data provides a more subtle example of complementarity. Despite the 
fact that the MAP error ellipse lies wholely within the SDSS error ellipse in Fig. the addition of SDSS provides 
substantially smaller error bars. As discussed in the combination of the CMB and galaxy power spectrum 
information yields a precise measurement of the Hubble constant h and flm independently of any low redshift GDM 
effects. The reason is that the physical extent of the sound horizon at recombination can be precisely calibrated 
from measurement of ftbh^ and flmh^ through the acoustic peak morphology (see Fig. || and [^). Measurement of 
this scale in redshift space isolates the Hubble constant; the aforementioned measurement of flmh^ in the CMB then 
returns Vim- Under the assumption of a flat universe, fig = 1 — flm is also well determined (see Fig. ^) and the 



^Note that in cases where strong model degeneracies and high experimental sensitivity coexist, 2000 modes are often required. 
For this reason, we do not quote errors for the Planck experiment alone; once combined with SN or galaxy survey data, its 
degeneracies are broken well enough that these small numerical errors are irrelevant. 

^A 68% confidence region for a two-dimensional ellipse extends to 1.52cr along each axis. We use this in all cases, although it 
is an overestimate in cases when the ellipse extends into unphysical or implausible regions of parameter space. 
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FIG. 5. Breaking the equation-of-state density degeneracy. The degeneracy exposed in Fig. ^ can be broken by combining 
CMB information (here from MAP) with SN or galaxy survey (here from SDSS) information. Plotted here are the 68% CL for 
the various experiments and combinations. Comparison of the two combinations leads to a sharp consistency test. Two fiducial 
models are shown: (a) Wg — —1/6 and (b) Wg — —1. While the CMB alone does well at Wg — —1/6, its degeneracy worsens 
considerably as Wg decreases toward —1. Note, however, that the complementary nature of the data sets and the ability to 
make consistency checks occurs in both cases. 



Cff = 0.03 


MAP(P)+SDSS 


MAP(P)+SN 


Planck(P)+SDSS 


Planck(P)+SN 


Wg 


a{wg) 


ai^lg) 


a{wg) 


a{ng) 


a{wg) 


a{Qg) 


a(wg) 


aiQg) 


-1/6 


0.015 (0.033) 


0.017 (0.028) 


0.024 


0.030 


0.009 (0.014) 


0.007 (0.010) 


0.014 


0.009 


-1/3 


0.027 (0.056) 


0.013 (0.022) 


0.028 


0.031 


0.016 (0.029) 


0.010 (0.019) 


0.020 


0.013 


-1/2 


0.047 (0.088) 


0.013 (0.022) 


0.041 


0.034 


0.022 (0.041) 


0.010 (0.019) 


0.021 


0.011 


-2/3 


0.074 (0.129) 


0.013 (0.022) 


0.063 


0.037 


0.029 (0.052) 


0.010 (0.020) 


0.023 


0.010 


-5/6 


0.108 (0.183) 


0.013 (0.022) 


0.091 


0.040 


0.037 (0.064) 


0.010 (0.019) 


0.026 


0.009 


-1 


0.126 (0.201) 


0.011 (0.018) 


0.125 


0.042 


0.033 (0.050) 


0.008 (0.016) 


0.027 


0.010 



TABLE I. Errors on Wg and fig upon combining data sets. SDSS assumes information out to fcmax ~ 0.2/iMpc ^ or 
O.l/iMpc"^ (parentheses). For Wg < —1/6, errors are insensitive to the sound speed of the fiducial model; here we use 
Coff = Vom « 0.2. 



angular diameter distance dAi^mh^ ,^g,Wg) may be used to extract Wg. Hence, despite the fact that galaxy surveys 
cannot determine these parameters by themselves, they can break the angular diameter distance degeneracy of the 
CMB and thereby allow the CMB to measure Wg and fig. The subtle nature of the degeneracy breaking requires a 
full joint analysis to uncover (c.f. |^] who reached more pessimistic conclusions from a separate examination of each 
data set). 

If the measurements pass the consistency test, we can combine all three sets of data. Even assuming no polarization 
information from MAP, the result is that the errors on Wg in the worst case of Wg = —1 become <j{wg) = 0.056, allowing 
for a sharp test for the presence of a cosmological constant. Failure to achieve consistency would indicate that one of 
our assumptions is wrong, e.g. Wg varies strongly with time or spatial curvature does not vanish fig 7^ 1 — flm- 

Finally, note that by marginalizing CcS, our results treat the clustering properties of the GDM as unknown and are 
thus conservative in the context of scalar-field quintessence models For example if Ccff is held fixed, the limits 
on Wg for the Wg — —1/6 model improve by ~ 30% for MAP-I-SDSS with or without polarization; gains are negligible 
near Wg — ~1 since Ccff has little observable effect there. 

In summary, for any value of Wg , the combination of CMB data with SN distance measures or galaxy surveys will 
provide reasonably precise measures of Wg and fig even considering the unknown clustering properties of the GDM. 
The comparison of these two combinations provides a sharp consistency test. 
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FIG. 6. Measurement of h and the equation of state degeneracy. Acoustic features in both the CMB (k^Pj) and matter/galaxy 
power (Pm) spectra are frozen in at last scattering. Once the CMB acoustic peaks are calibrated in real space from the peak 
morphology measurement of Qmh^ and ftbh^, sliding the galaxy power spectrum in redshift space {h Mpc~^) until the features 
"match" determines h. This test is unaffected by and therefore insensitive to late-time dynamics from the GDM component. 
However, once h is determined, Qg — 1 ~ Qrn follows from the CMB measurement of flmh^- The angular diameter distance 
measurement from the CMB then determines w„. 




FIG. 7. Significance of bounds on c^s as a function of the equation of state Wg of the fiducial model. Plotted is the formal 
significance with which the simplest scalar-field models (coh = 1) may be excluded if the true model is c^g — \/0.03 ~ 0.2. The 
formal errors may overestimate the significance as discussed in the text. The power to distinguish sound speed effects decreases 
sharply as one approaches the cosmological constant case of uig — —1. 
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no priors: MAP+SDSS MAP(P)+SDSS Planck+SDSS Planck(P)+SDSS 

0.03 0.13 (0.13) 0.06 (0.06) 0.12 (0.12) 0.04 (0.04) 

0.1 0.46 (0.47) 0.14 (0.15) 0.42 (0.43) 0.10 (0.10) 

0.3 1.52 (1.55) 0.54 (0.54) 1.47 (1.47) 0.37 (0.37) 
cr(ln/3) = 0.1: 

0.03 0.06 (0.06) 0.06 (0.06) 0.06 (0.06) 0.04 (0.04) 

0.1 0.18 (0.20) 0.14 (0.15) 0.18 (0.18) 0.10 (0.10) 

0.3 0.75 (0.83) 0.52 (0.53) 0.72 (0.72) 0.37 (0.38) 



TABLE II. Errors on c^g as a function of its fiducial value for Wg = —1/3. SDSS assumes fcmax = 0.2/i Mpc~^ {O.lh Mpc~^). 
P denotes the inclusion of CMB polarization information. Upper: no priors on P; lower: prior of a(ln/3) = 0.1. The prior helps 
isolate r if polarization information is unavailable and hence brings the limits with and without polarization closer. 



V. CONSTRAINING THE SOUND SPEED 



If the tests of the last section determine that the equation of state of the exotic component is in the range — 1 < 
Wg < 0, we will have discovered a new form of matter. It then becomes interesting to explore its properties in order 
to search for a suitable particle physics candidate. The simplest candidate is a slowly-rolling scalar field, also known 
as "quintessence" [p^l^^Tsf . The hallmark of such a candidate is that its effective sound speed is simply the speed of 
light (coff = 1), i.e. it is a maximally stable form of matter. It furthermore has Cvis — 0. Can the clustering properties 
of the GDM be measured well enough to distinguish a scalar field component from alternate candidates for the exotic 
matter? 

Stabilization of perturbations may occur through a finite effective sound speed CcS as it does for a real scalar field 
or through viscosity Cvis as some defect-dominated models suggest ||l^,|l8|. Because the two are largely degenerate, 
we take Coff as a proxy for both - one actually determines the combination of Coff and Cvis that fixes the GDM sound 
horizon. 

The sound speed CcS will only be well constrained if the GDM sound horizon at Zg is sufficiently small that 
cancellation of the ISW effect varies strongly with Ccs (see Fig. ||) or if the features in the matter power spectrum lie 
on scales accessible to galaxy surveys. Both of these considerations favor fiducial models with low Coff. SN distance 
measures have no dependence on . In Table || (upper), we show the errors on Ccff as its fiducial value increases in 
an Wg = —1/3, fig = 0.65 fiducial model. As usual, both these and other parameters are marginalized when quoting 
errors on c^g. As expected, the error a{c^g) increases sharply as Cofr — > 1. 

The strong variation of cr(c^jf) with c^g itself makes it difficult to estimate the significance at which two models 
can be separated. For example, if we were to take cr(c^g = 0.3) to infer that it is distinguishable from c^g = at 
only 0.3/1.52 = Q.2a from MAP+SDSS, we would be incorrect since a model with a smaller value c^ff = 0.03 is 
distinguishable from zero at a higher level. Conversely, the ability to distinguish a fiducial model with Ccff = from 
one with CoS > is always overestimated. This problem reflects the limitations of the Fisher matrix technique caused 
by its infinitesimal expansion of the likelihood function. 

To address this issue, we take an intermediate value of Coff = y/0.03 ~ 0.2 and ask how well we can reject the scalar 
field hypothesis of Ccs — 1. We plot in Fig. |^, the number of standard deviations by which the true sound speed is 
separated from the scalar field value ["cr's" = (1 — c^ff)/CT(c^ff)]. Although this formal significance still overestimates 
the true significance, the qualitative result is clear. If Wg > —1/2, CMB and galaxy survey data will be able to place 
interesting constraints on the sound speed. As Wg decreases to —1, the effects of clustering in the GDM vanish leaving 
no significant constraint on Ccff. A more complete exploration of the likelihood function would yield more precise 
limits but is computationally time consuming; we defer such an analysis until Wg is measured and shown to be in this 
range. 

Fig. implies that good limits on Ccff depend on polarization information. This is because reionization effects are 
nearly degenerate with ISW effects from Coff. We show this in Fig. ||. Polarization information isolates r from the 
feature generated by Thomson scattering of anisotropic radiation present at large scales during reionization. However, 
because of foregrounds and systematics likely in the large-angle polarization data, it is interesting to see whether any 
other information can break this degeneracy. The main effect of r is to reduce the small-angle anisotropics in the 
CMB uniformly. If the intrinsic amplitude can be calibrated by the galaxy survey data, r could be measured. With 
the growth function and other transfer function effects under control, the remaining obstacle is the unknown bias 
factor b. This can be measured on large scales through redshift-space distortions. Since flm is well-constrained by 
the combination of CMB and galaxy survey data (§ IV), the constraint on f3 — /h from these distortions supply 



10 




Ceff2 



FIG. 8. The T — c^g degeneracy, r and c^g are nearly degenerate given only temperature information. With either polarization 
information on the CMB side [MAP(P)] or a prior of (j(ln/3) — 0.1 from redshift space distortions on the galaxy power spectrum 
side [SDSS(/3)], the de generacy is broken allowing better isolation of c^q. 



information on the bias. Taking a conservative prior of (T(ln/3) — 0.1 |36|, the normalization determination breaks the 
Cgfj — T degeneracy almost as effectively as polarization information. We quantify this in Tab. ^ (lower), where the 
errors on c^g with and without polarization are made more comparable with this conservative prior on 

In summary, interesting constraints on the clustering properties of the exotic component will be available if Wg ^ 
— 1/2 as long as we have either polarization data from the CMB or redshift-space distortion information from galaxy 
surveys. 



VI. DETECTING THE NEUTRINO BACKGROUND RADIATION 



If the equation of state of the GDM is determined to be = — 1, then the only possibility is a cosmological 
constant. In this case, the basic aspects of structure formation are so simple that subtle effects in the dark sector can 
be uncovered. The remaining dark matter in the universe is the neutrino background radiation (NBR). How well can 
we detect its presence? The issue is somewhat more subtle than it initially appears due to a degeneracy in the CMB 
acoustic peaks. We shall see that detecting the neutrino background radiation requires detecting its fluctuations, in 
particular its anisotropics. 



A. Matter-Radiation Degeneracy 



Given that CMB anisotropics are generally sensitive to changes in the expansion rate at high redshift, one might 
think the radiation content of the universe could be measured precisely. Indeed the matter density flmh^ can be 
measured to a{il.„ih^) = 0.02 by the MAP satellite (without polarization) if the radiation is taken to be fixed. The 
problem is that what the CMB best measures is the matter-radiation ratio, not the matter or radiation density 
individually. 

The GDM parameterization can be used to explore this degeneracy and more generally deconstruct the information 
contained about the NBR in the CMB. We know that the matter-radiation degeneracy arises because of the way the 
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FIG. 9. Matter-Radiation degeneracy. Plotted here are the 68% CL for the various experiments, (a) Ignoring the information 
supplied by the anisotropy of the NBR by marginalizing over Cvis, we find a strong degeneracy between the neutrino number 
Ni, and the matter density Qmh^ along a line of constant matter-radiation ratio pm/pr in both the CMB and galaxy survey 
data. Combining CMB information from MAP and galaxy information from SDSS breaks the degeneracy somewhat, (b) 
Utilizing the information in the anisotropies by assuming Cvis = y^1/3 shrinks and rotates the error ellipses leading to better 
complementarity and tighter errors. 



background expansion rate scales with these parameters [see Eq. (^]. Fluctuations in the matter and radiation break 
this degeneracy. We can use the GDM parameterization to separate the information on the background expansion 
rate from that of the fluctuation properties. 

As shown in and H, the NBR is accurately modeled by a GDM component with Wg = c^g — c^jg = 1/3 and 
density Qg — 5.63 x 10^"h^'^N^. A^i, = 3 in the fiducial model. Thus N^, and Wg determine the background properties 
whereas Ceff and Cvis control the fluctuation properties, with Cvis controlling the anisotropic stress of the NBR. The 
anisotropic stress is proportional to the quadrupole anisotropy of the NBR. Thus Cvis = (with CoS — ^1/3) represents 
a component with the same background properties as the NBR but with no anisotropies. 

How does ignorance of the properties of the NBR affect the determination of the matter density Om/i^? If we 
allow Nu and Cvis to vary so as to eliminate the information provided by the density and anisotropy of the radiation 
component, the error ellipses of Fig. ^ reveal a matter-radiation degeneracy, i.e. they are elongated along the line of 
constant matter-radiation density ratio. The degeneracy affects both CMB measurements and galaxy surveys alike. 
Here and in the remainder of this section, we keep Wg and Ceff fixed while varying 17^, Cvis, and the other cosniological 
parameters including a cosmological constant ^Ijy. The MAP errors on ^Imh'^ are degraded from cr(ilm/i^) = 0.02 
to 0.16 when is allowed to vary. The baryon-to-photon ratio flbh'^ remains well- measured. Some leverage in a 
flat universe is provided by the fact that the actual matter density comes into the angular diameter distance for the 
CMB. Indeed, if one fixes the other parameters that go into the angular diameter distance, in this context JIa, then 
the CMB does place tight constraints on flmh^ and N^, separately However, in the general case, the angular 

diameter distance degeneracy prevents a precise measurement of flmh^ by these means. 

Combining CMB anisotropies and galaxy power spectrum information, which both suffer from the matter-radiation 
degeneracy individually, restores tight error bars on ^mh"^ even using only temperature information from MAP (see 
Tab. III). Here the additional information on rih/rim from baryonic features in the galaxy power spectrum along 
with the precise measurement of r2f,/i^ from the CMB constrains Q.mh? ■ This is another example in which the 
complementary nature of the CMB and galaxy survey data helps in a subtle way. 
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MAP+SDSS MAP(P)+SDSS Planck+SDSS Planck(P)+SDSS 

assumption a{flrn,h^) o-{Nv) a{Qrn.h^) o'{Nv) aiyimh^) u{N,y) (j{Q,mh?) ct{Nv) 

unknown cjj, 0.026 (0.046) 1.24 (2.26) 0.023 (0.036) 1.12 (1.92) 0.006 (0.006) 0.30 (0.43) 0.004 (0.004) 0.21 (0.23) 

fixed c^i^ 0.007 (0.024) 0.44 (1.59) 0.006 (0.022) 0.43 (1.44) 0.003 (0.005) 0.23 (0.43) 0.003 (0.003) 0.17 (0.20) 

TABLE III. Errors on Q,mh^ and N^, with and without information from NBR anisotropies. SDSS assumes fcmax = 0.2/i 
Mpc~^ (fcmax = O.l/i Mpc~^). P denotes the inclusion of CMB polarization information. 



B. Limiting the Number of Neutrinos 



In the standard scenario, the fluctuations of the NBR are not unknown; they are fixed through the properties of the 
neutrinos and gravitational instability. These fluctuations in the NBR further breaks the matter-radiation degeneracy. 
Like the CMB itself, the NBR carries temperature anisotropies (see |2^] for the full angular power spectrum). In 
particular, the quadrupole anisotropy of neutrinos alters the gravitational potentials that drive acoustic oscillations 
(see Fig. I). 

By fixing Cvis — ^^1/3, we can ask how well the number of neutrino species in the NBR can be measured 
under the standard assumptions. This is appropriate for either flavored or sterile neutrinos if their mass is sufficiently 



small {nil, ^ 0.1 eV). We show the results in Fig. ^b and Tab. III. The CMB and galaxy survey ellipses in the 
Njj-^rn.h?' plane shrink and rotate in opposite senses from the case of a marginalized Cvis (Fig. ||a). This enhances 
the complementary nature of CMB and galaxy survey data giving a substantial improvement when the data sets are 
combined if only MAP data is available. 

As the errors in Tab. HI are comparable to those achievable from big bang nucleosynthesis (BBN) |4^], CMB 
and galaxy surveys should provide a powerful consistency check on BBN and a constraint on additional neutrino 
species populated around recombination. Unfortunately, testing percent-level differences in the NBR temperature 
(represented here by a change 5N,j/N,y — A6Ti,/Ti,, see Eq. Q) due to the details of their decouphng p7|-p9| seems 
out of reach even if we combine all of our precision tests. 



C. Detecting Neutrino Anisotropies 



Anisotropies in the NBR are predicted by the gravitational instability paradigm and are potentially observable 
through their effect on CMB anisotropies. In the GDM model for the NBR, the anisotropies are determined by the 
viscosity parameter Cvisi constraints on this parameter tell us how well anisotropies in the NBR can be detected. The 
role of the viscosity parameter Cvis in breaking degeneracies in the last section suggests that the anisotropies in the 
neutrino background radiation may themselves be detectable. Unfortunately, with the CMB alone, the effect is strongly 
degenerate with those of other parameters. Even though we fix the other properties of the GDM {wg = c^g = 1/3, 
Ni, = 3), changes in Cvis can be mimicked by changes in the normalization and tilt of the spectrum at small angles. 

By adding in galaxy survey data, c^jg = 1/3 (NBR anisotropies) and c^jg = (no anisotropies) are separated by 
3.5a from MAP-I-SDSS. The significance improves to 8.7a with Planck. 

How sensitive is the measurement to the underlying assumptions about the data set and model space? The loss of 
polarization information does not significantly affect these limits. On the other hand, if we take the more conservative 
fcmax = 0.1/iMpc^^ for the galaxy surveys, the significance decreases to 1.2a {7.3a) for MAP-KSDSS (Planck -l-SDSS). 
Perhaps more important, the MAP-I-SDSS result does depend on prior knowledge that N^, w 3. Fortunately, even 
assuming only very conservative constraints of a{N,^) = 1.0 from big bang nucleosynthesis, allows a detection at 
2.0a {7.1a) for MAP-I-SDSS (Planck-|-SDSS). These results imply that NBR anisotropies can be detected with high 
significance at least by the Planck satellite, even under conservative assumptions. 



VII. CONCLUSIONS 



With the wealth of precision cosmological measures that are becoming available, we should soon be in the position 
to identify all of the cosmologically important components of the universe - including any dark components that may 
be present. Not only do CMB anisotropies, high-redshift objects, and galaxy surveys probe different aspects of the 
cosmology, but they can work together to uncover a complete and consistent cosmological picture. 
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FIG. 10. Detecting anisotropies in the neutrino background radiation. Anisotropics in the NBR are detectable if a model 
with c^ij. = 1/3 can be distinguished from one with cJig = 0. Degeneracy with other parameters such as the tilt n prevent the 
CMB information provided by MAP from detecting the anisotropies. Adding galaxy survey information from SDSS breaks the 
degeneracy and allows a statistically significant detection. We have assumed here that N,j is fixed at 3. 



We have shown here how the combined power of these data sets can determine the properties of the dark components. 
Assuming that the prehminary indications from SN data that the missing component is not merely spatial curvature 
are confirmed, the first step will be to determine the equation of state of the exotic component. The task is non-trivial 
due to a degeneracy with its density in determining the expansion rate. The degeneracy is broken by combining the 
CMB data with SN distance measures, galaxy surveys, or any other measurement that can constrain = 1 — or 
h. By further combining any of these pairs of data sets, we create powerful consistency tests. Note that these tests 
work even near a cosmological constant model Pg/ Pg = Wg = —1, where the degeneracy in the CMB alone is at its 
worst. 

Should the equation-of-state measurement rule out a cosmological constant, wc will want to study the other prop- 
erties that identify this exotic component. Its clustering scale is accessible in the CMB and galaxy surveys as features 
in their power spectra at large scales. We have parameterized this with a sound speed and have shown that as long 
as the equation of state is sufficiently different from a cosmological constant {wg ^ —1/2), one can distinguish a 
maximally stable component {ccS = 1) from a component with CeS ^ 0.2. Distinguishability increases substantially as 
Pg/pg increases, reaching CcS ^ 0.6 at = — 1/6. Such limits are interesting since the simplest physically motivated 
exotic component, a slowly-rolling scalar- field quintessence, is maximally stable with Coff = 1. 

We have only considered models where the equation of state varies in time sufficiently slowly to be replaced by 
some suitably averaged but constant Wg. If the GDM sector involves stronger temporal variation, to what extent will 
upcoming data sets be able to constrain the possibilities? As noted in § [V , acoustic features in the CMB anisotropies 
when combined with those in the galaxy power spectrum can be combined to measure flm- This measurement requires 
that the dark components at high redshift such as the NBR are known but makes no assumptions whatsoever about 
the low- redshift behavior of Wg. With the present-day value of fig = 1 — r2,„ known to fair accuracy, there are a 
number of observational handles on Wg as a function of time. The location of the CMB acoustic peaks determines 
the angular diameter distance to high redshift. Mid-redshift supernovae constrain the luminosity distance to z w 0.5, 
although with very large samples one may even extract some redshift dependence. Measurements of the normalization 
of the power spectrum on scales below the GDM sound horizon will constrain the GDM-modified growth rate. The 
z = normalization may be measured from abundances of rich clusters and from the galaxy power spectrum given 
a measurement of galaxy bias from redshift distortions or peculiar velocity data sets. The normalization at higher 
redshift can be estimated from the statistics of the Lyman a forest, damped Lyman a systems, and high-redshift 
clusters. Hence, although the most general equation of state is described by a free function of redshift, there are 
actually a number of robust observational handles on its behavior! 
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A similar analysis shows that even if the universe contains both an accelerating component and a non-vanishing 
spatial curvature, which in many respects resembles a universe with Wg varying from —1/3 to —1, the combination of 
information from different redshifts will give us leverage on the two separately. The situation is actually even more 
favorable since the geometrical aspects of curvature enter strongly into the angular diameter distance as measured by 
the CMB. 

Should the measurement of the equation of state confirm the relative simplicity of a cosmological constant, we 
will be able to probe in detail the remaining dark component, the neutrino background radiation. Detection of 
the neutrino background radiation through the CMB suffers from the fact that a change its energy density may be 
compensated with a change in the matter density up to effects due to the presence of fluctuations. We have shown 
that by combining CMB and galaxy survey data, this degeneracy can be broken and yields constraints competitive 
with those from big bang nucleosynthesis. Furthermore, the combination of CMB and galaxy survey data should 
provide the first detection of anisotropics in the neutrino background radiation. The detection of these anisotropics, 
predicted by the gravitational instability paradigm for structure formation, would represent a triumph for cosmology. 
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